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Abstract 

We discuss the problem of ultra high energy particles propagation in astrophysical backgrounds. We present two difTerent 
computational schemes based on both kinetic and Monte Carlo approaches. The kinetic approach is an analytical 
computation scheme based on the hypothesis of continuos energy losses while the Monte Carlo scheme takes into account 
also the stochastic nature of particle interactions. These schemes, that give quite reliable results, enable the computation 
of fluxes keeping track of the different primary and secondary components, providing a fast and useful workbench to 
study Ultra High Energy Cosmic Rays. 
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Introduction and 
Conclusions 



Ultra High Energy Cosmic Rays (UHECR) are the 
niost energetic particles observed in nature, with en- 
ergies up to few X 10^° eV. Nowadays the experimen- 
tal study of UHECR is conducted by three different 
experiments: Auger in Argentina, HiRes and Tele- 
scope Array in the USA. 

The propagation of UHECR from the source to 
the observer is conditioned by their interactions with 
astrophysical backgrounds: the Cosmic Microwave 
Background (CMB) and the Extragalactic Back- 
ground Light (EBL). Understanding the key features 
of the propagation is of paramount importance to 
interpret experimental observations paving the way 
for the discovery of the astrophysical origin of these 
fascinating particles. 

Several features of the observed spectrum can 
be directly linked to the chemical composition of 
UHECR and to their sources (Greisen 1966, Zatsepin 
& Kuzmin 1966, Berezinsky et al. 2006a, Aloisio et. 
al. 2007, Aloisio & Boncioli 2011). Among such 
features particularly important is the Greisin, Zat- 
sepin and Kuzmin (GZK) suppression of the flux, an 
abrupt depletion of the observed proton spectrum, 
arising at energies £' ~ 5 x 10^^ eV, due to the in- 
teraction of UHE protons with the CMB radiation 
field (Greisen 1966, Zatsepin & Kuzmin 1966). The 
GZK suppression, as follows from the original pa- 
pers, is referred to protons and it is due to the photo- 
pion production process on the CMB radiation field 



(p+jCAiB — > TT+p). In the case of nuclei the expected 
flux also shows a suppression at the highest energies 
that, depending on the nuclei specie, is due to the 
photo-disintegration process on the CMB and EBL 
radiation fields {A -\- ')cmb,ebl {A — nN) A- nN) 
(Aloisio et. al. 2012a). Another important fea- 
ture in the spectrum that can be directly linked with 
the nature of the primary particles and their origin 
(galactic/extra-galactic) is the pair-production dip 
(Berezinsky et al. 2006a, Aloisio et al. 2007). This 
feature is present only in the spectrum of UHE extra- 
galactic protons and, as the GZK, is a direct conse- 
quence of the proton interaction with the CMB radia- 
tion field, in particular the dip brings a direct imprint 
of the pair production process p-|-7cA/B — >■ p+e^-\-e~ 
suffered by protons. 

From the experimental point of view the situ- 
ation is far from being clear with different experi- 
ments claiming contradictory results. The HiRes ex- 
periment, not anymore taking data, showed a proton 
dominated spectrum till the highest energies (HiRes 
collaboration 2008, 2010) while the Auger observa- 
tions show an heavy mass composition at energies 
> 4 X 10^8 eV (Auger collaboration 2010). 

This puzzling situation, with different experi- 
ments favoring different scenarios, shows once more 
the importance of a systematic study of UHECR 
propagation in astrophysical backgrounds. In the 
present paper we will review the main points of two 
alternative computation schemes which enable the 
determination of the fiuxes expected on earth fix- 
ing the injection spectrum and the distribution of 
sources. These two schemes are based on different 
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approaches to model particle interactions with back- 
grounds: the continuum energy losses (CEL) approx- 
imation, which is the base for the kinetic approach, 
and the Monte Carlo (MC) technique. 

As we will discuss in the following these two dif- 
ferent schemes give reliable results that, in the frame- 
work of the different assumptions, agree each other 
offering a suitable theoretical framework to study ex- 
perimental results unveiling the intimate nature of 
UHECR. 



2 Kinetic Equations 

The main assumption under which the kinetic the- 
ory is build is the CEL approximation (Berezinsky 
et al 1990), through which particle interactions are 
treated as a continuum process that continuously de- 
pletes the particles energy. 

UHECR propagating through astrophysical back- 
grounds suffer different interaction processes 

• protons - UHE protons interact only with the 
CMB radiation field giving rise to the two pro- 
cesses of pair production and photo-pion pro- 
duction. Both of these reactions can be treated 
in the CEL hypothesis. 

• nuclei - UHE nuclei interact with the CMB and 
EBL radiation fields, suffering the process of 
pair production, for which only CMB is rel- 
evant, and photo-disintegration, that involves 
both backgrounds CMB and EBL. While the 
first process can be treated in the CEL hypoth- 
esis, being the nuclei specie conserved, the sec- 
ond cannot be, producing a change in the nu- 
cleus specie. Following Aloisio et. al. 2012a, in 
the framework of the kinetic approach, we will 
treat the photo-disintegration process as a " de- 
caying" process that simply depletes the flux of 
the propagating nucleus. 

Taking into account all energy losses processes we 
can describe the propagation of protons and nuclei 
through kinetic equations of the type: 



dnp{T, t) d 



dt 



dT 



[6p(r,0np(r,t)] -Qp(r,o (i) 



QA{T,t) 

where n is the equilibrium distribution of particles, b 
are the energy losses (adiabatic expansion of the Uni- 
verse and pair/photo-pion production for protons or 
only pair-production for nuclei) Q is the injection of 
freshly accelerated particles and, in the case of nuclei, 
also the injection of secondary particles produced by 
photo-disintegration (see below). 



The energy losses b for protons or nuclei depend 
only on the CMB field and in the CEL hypothesis can 
be computed analytically (Berezinsky et. al. 2006a, 
Aloisio et. al . 2007, Aloisio et. al. 2012a). 

The second process that affects nuclei propaga- 
tion is photo-disintegration over CMB and EBL back- 
grounds. This process is treated as a decaying pro- 
cess that depletes the fiux of nuclei, it enters in 
the kinetic equation (see equation ([2])) through a 
sort of "life-time" of the nucleus under the photo- 
disintegration process. This "life-time" corresponds 
to the mean time needed to a nucleus of Lorentz fac- 
tor r and atomic mass number A to lose, at least, 
one of its nuclcons: 
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where (7{er,A) is the photo-disintegration cross- 
section and i^(er) is the molteplicity associated with 
this process, namely the average number of nucleons 
extracted from the nucleus by a single interaction and 
n-bkg = n-cMB + nEBL- The dependence on red-shift 
of TA directly follows from the evolution with red- 
shift of the background photon densities ucmb and 
nEBL- In the case of CMB this dependence is known 
analytically while for the EBL one should refer to 
evolution models (in our computations we have used 
the model by Stecker et. al. 2006). 

One important feature of the photo-disintegration 
process is that it starts to contribute to the propa- 
gation of nuclei at a Lorentz factor that is almost in- 
dependent of the nuclei specie Tcr — 2 x 10^ (Aloisio 
et al. 2012a). This is an important general char- 
acteristic of nuclei photo-disintegration process from 
which we can immediately deduce the dependence 
on the nucleus specie of the energy corresponding 
to the photo-disintegration suppression of the flux: 
-E^j = AmN^cr being A the atomic mass number of 
the nucleus and ruN the proton mass. From this ex- 
pression for it is evident how the flux behavior 
could bring informations on the chemical composi- 
tion of the UHECR, in the case of Helium {A — 4) 
the suppression is expected around energies E ~ 10^^ 
eV while in the case of Iron {A = 56) the suppression 
is expected at higher energies E ~ 10^" eV. 

Let us discuss now the generation function 
QA{^,t) in the right hand side of Eq. ([2]). One 
should distinguish among primary nuclei, i.e. nuclei 
accelerated at the source and injected in the inter- 
galactic space, and secondary nuclei and nucleons, 
i.e. particles produced as secondaries in the photo- 
disintegration chain. In the case of primaries the in- 
jection function is an assumption of the source model, 
while the injection of secondaries should be mod- 
eled taking into account the characteristics of the 
photo-disintegration process. The dominant process 
of photo-disintegration is the one nucleon (TV) emis- 
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sion, namely the process {A + 1) + ^bkg A + N, 
this follows directly from the behavior of the photo- 
disintegration cross-section (see Aloisio et al. 2012a 
and references therein) that shows the giant dipole 
resonance corresponding to one nucleon emission. 
Moreover, at the typical energies of UHECR {E > 
10^^ eV) one can safely neglect the nucleus recoil so 
that photo-disintegration will conserve the Lorentz 
factor of the particles. Therefore the production rate 
of secondary A— nucleus and A— associated nucleons 
will be given by 



QAiT,z) 



nA+i{T,z) 
TA+i{r,z) 



(4) 



where ta+i is the photo-disintegration life-time of the 
nucleus father (^+1) and ua+i is its equilibrium dis- 
tribution, solution of the kinetic equation ©. 

Using equation ^ we can build a system of cou- 
pled differential equations that starting from pri- 
mary injected nuclei (^o) follows the complete photo- 
disintegration chain for all secondary nuclei {A < Aq) 
and nucleons. Clearly secondary proton^ propaga- 
tion will be described by the proper kinetic equation 
([1]) with an injection term given by Q. The solution 
of the kinetic equation for protons and nuclei can be 
worked out analytically. In the case of protons: 



dz' 



{l + z')H{z') 



(5) 



being Qp the injection of primary protons or sec- 
ondary protons (equation @) and T' — r'(r,z) is 
the characteristic function of the kinetic equation 
(Aloisio et al. 2012a). In the case of nuclei: 



nA(r, z) 



(1 



z')i?(^') ' dV 

(6) 

being, again, Qa the injection of primary or sec- 
ondary ([4]) nuclei. The exponential term in Eq. (jG)) 
represents the survival probability during the propa- 
gation time t' — t for a nucleus with fixed A and can 
be computed according to Aloisio et al. 2012a. The 
derivative term dT' /dT present in both solutions ^ 
and dH is analytically given (Aloisio et al. 2012a). 



3 Monte Carlo 

The kinetic approach outlined above neglects inter- 
actions fluctuations considering an (average) contin- 
uum loss of energy suffered by particles. This approx- 
imation in the case of protons, has a limited effect 
on the flux computation only at the highest energies 
{E> 100 EeV) (Berezinsky et al. 2006a, Berezinsky 
et al. 2006b, Aloisio et al. 2007). 
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Figure 1: Flux of iron and secondary nuclei (A=50, 40, 
30, 20, 10) at z = in the case of pure iron injection at 
the source with a power law injection index 7 = 2.2. Full 
squares correspond to the SimProp result (Aloisio et al. 
2012b) while continuous lines correspond to the solution 
of the nuclei kinetic equation of Aloisio et al. 2012a. 

In order to evaluate the effects of fluctuations on 
the expected nuclei flux we have build a computation 
scheme alternative to the kinetic one, that uses the 
MC technique to simulate nuclei interactions. First 
of all, let us remark that fluctuations could be rele- 
vant only in the case of nuclei photo-disintegration, 
this follows from the fact that the pair-production 
process involving nuclei can be considered as an inter- 
action process of the inside nucleon, therefore being 
fluctuations in the protons pair-production irrelevant 
(Berezinsky et al. 2006b) the same holds for nuclei. 
The MC simulation scheme we have developed Sim- 
Prop (Aloisio et al. 2012b) is mono-dimensional: it 
does not take into account spatial distributions tag- 
ging sources only through their distance from the 
observer (red-shift). The MC simulation propagates 
particles in steps of red-shift following the injected 
nucleus, secondary nuclei and protons produced at 
each photo-disintegration interaction and calculating 
their losses up to the observer, placed at red shift 
zero. The nuclear model on which SimProp is based 
is the same used for the kinetic approach (see Aloisio 
et al. 2012a, 2012b and references therein). The 
stochastic nature of the nuclei photo-disintegration 
process is modeled through the survival probability 
of a nucleus of atomic mass number A and Lorentz 
factor F 



P(F, z) = cxp 



1 



ta(X,z') 



dt 
dz' 



dz 



(7) 



where z and z* are the values of the redshift of the 
current step (from z* to z). 

The SimProp code is designed in such a way that 
any red-shift distribution of sources and any injection 
spectrum can be simulated. This is achieved draw- 
ing events from a flat distribution in the red-shift 



^Neutrons decay very fast into protons, so we will always refer to secondary protons. 
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of the sources and of the logarithm of the injection 
energy. Once the event is recorded at z = the ac- 
tual source/energy distribution is recovered through 
a proper weight attributed to the event (Aloisio et 
al. 2012b). We will compare now the spectra ob- 
tained using SimProp (Aloisio et al. 2012b) with 
those calculated solving the kinetic equation asso- 
ciated to the propagation of nuclei (Aloisio et al. 
2012a). To pursue such comparison, a pure iron in- 
jection with a power law injection of the type oc E~'^ 
with 7 = 2.2 have been assumed. The sources have 
been assumed to be homogeneously distributed in the 
red-shift range < z < 3. In figure [1] the fluxes ex- 
pected at z = are shown for iron and secondary 
nuclei produced in the photo-disintegration chain suf- 
fered by primary injected irons. The points refer to 
the SimProp results while the continuous lines to the 
fluxes computed in the kinetic approach. A good 
agreement between the two schemes is clearly visible 
in figure [T] At the highest energies the path-length 
of iron nuclei is very short. Therefore, to achieve a 
good sampling in the MC simulation, higher statis- 
tics is needed; this is the reason for larger errors bars 
in the SimProp results at the highest energies. Let 
us conclude discussing why it is useful to go beyond 
the kinetic approach. The kinetic approach has the 
important feature of being analytical: fluxes are com- 
puted mathematically solving a first principles equa- 
tion (Aloisio et al. 2012a). This means that the fiux 
of primaries and secondaries is expressed in terms of 
several integrals that can be computed numerically, 
once the injection spectrum and the sources distri- 
bution are specified. In particular, the flux of sec- 
ondary nuclei and nucleons produced in the photo- 
disintegration chain of a primary Aq is obtained by 
the numerical computation of Aq nested integrals and 
this computation should be repeated each time the 
hypothesis on sources (injection and distribution) are 
changed. This computation, while it is always feasi- 
ble numerically, takes some time that can be substan- 
tially reduced using a MC computation scheme. This 
follows by the fact that, as discussed above, within 
the SimProp approach it is possible to simulate differ- 
ent source distributions and injection spectra without 
repeating the overall propagation of particles. In this 
sense the MC approach presented here, which is the 
minimal stochastic extension of the kinetic approach, 
provides a faster computation scheme. 
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DISCUSSION 

CARLO GUSTAVINO: The difference between 
Auger and the other experiments can be due to the 
fact they are looking from different hemispheres? 

ROBERTO ALOISIO This is an hypothesis that 
was recently put forward. I personally do not believe 
in such explanation because of the simple reason that 
at energies around 2-^3 x 10^^ eV, where already the 
difference between Auger and HiRes starts, the uni- 
verse visible in UHECR has a huge scale of the order 
of Gpc. Therefore it is very unlikely to have dif- 
ferences between observations carried out from the 
southern and northern hemispheres. 



